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1.  Molecules  studied  in  this  work:  1,  4-methyl-2-oxetanone;  2,  6,8-dioxabicyclo[3.2.1] 

octane;  3,  and  l,7,7-trimethylbicyclo[2.2.1]  heptan-2-one  .  1 


2.  Mid-ER  absorption  and  VCD  spectra  of  R  -  (+)  -  1.  Experimental  spectra  (a  and  c) 
from  Devlin  et  al.  (to  be  published).  Calculated  spectra  (b  and  d)  use  6-3 IG* 

B3LYP/DFT  force  field  and  APTs;  AATs  are  6-3 IG*  semi-B3LYP/DFT.  Calculated 
spectra  use  Lorentzian  band  shapes  with  y  =  4.0  cm~^  for  all  bands.  Note  that  some 
of  the  differences  between  calculated  and  experimental  absorption  and  VCD  intensities 
originate  in  deviations  in  experimental  band  widths  from  y  =  4.0  cm“^  The 

assignment  of  fundamentals  is  indicated;  see  also  Table  1 .  7 

3.  Mid-IR  absorption  and  VCD  spectra  of  (IR,  5S)  -  (+)  -  2.  Experimental  spectra 
(a  and  c)  from  Eggimann,  Shaw,  and  Wieser  (1991)  and  Eggimann  et  al.  (1993). 

Calculated  spectra  (b  and  d)  use  6-31G*  B3LYP/DFT  force  field  and  APTs;  AATs 

are  6-3  IG*  semi-B3LYP/DFT.  y  =  4.0  cm“^  for  all  bands.  The  assignment  of 

fundamentals  is  indicated;  see  also  Table  2 .  9 
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force  field  and  APTs;  AATs  are  6-31G*  semi-B3LYP/DFT.  y  =  4.0  cm~^  for  all 

bands.  The  assignment  of  fundamentals  is  indicated;  see  also  Table  3 .  11 

5.  Mid-IR  absorption  and  VCD  spectra  of  R  -  (+)  -  L  Spectra  a  and  d  are  obtained 
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elimination  of  all  nonfundamental  bands.  Calculated  spectra  (b,  c,  e,  f,  and  g)  use 
B3LYP/DFT  force  field  and  APTs;  AATs  are  semi-B3LYP/DFT.  In  b  and  e,  the 
basis  set  is  T^P;  in  c  and  g,  it  is  6-3  IG*.  In  f,  the  force  field  and  APTs  are 
6-3 IG*,  while  the  local  AATs  are  TZ/2P.  y  =  4.0  cm~^  for  all  bands  in  all 

calculated  spectra .  14 

6.  Mid-IR  absoiption  and  VCD  spectra  of  R  -  (+)  -  L  Spectra  a  and  d  are  the 
fundamental  experimental  spectra  from  Figures  5a  and  5d.  Calculated  spectra 

(b,  c,  e,  and  f)  use  the  6-3 IG*  B3LYP/DFT  force  field.  In  b  and  e,  the  APTs  are 

6-31G*  B3LYP/DFT  and  the  AATs  are  6-31G*  semi-B3LYP/DFT.  In  c  and  f.  the 

APTs  and  AATs  are  6-3 IG*  SCF.  y  =  4.0  cm”^  for  all  bands  in  all  calculated 

spectra  .  16 

7.  Mid-IR  absorption  and  VCD  spectra  of  R  -  (+)  -  1.  Spectra  a  and  e  are  the 
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1.  INTRODUCTION 


We  report  a  major  advance  in  the  ab  initio  theoretical  prediction  of  vibrational  circular  dichroism 
(VCD)  spectra.  This  advance  pennits  useful  predictions  of  VCD  spectra  for  much  larger  molecules  than 
heretofore  accessible.  As  a  result,  the  practical  utility  of  VCD  spectroscopy  in  smdying  the 
stereochemistry  of  chiral  molecules  is  greatly  enhanced. 


Specifically,  we  report  the  first  ab  initio  calculations  of  VCD  spectra  based  on  harmonic  force  fields 
obtained  using  density  fimctional  theory  (DFT).  DFT  is  increasingly  used  in  the  calculation  of  molecular 
properties  (Ziegler  1991;  Labanowski  and  Andzelm  1991).  Its  growing  popularity  is  due  to  (1)  the 
development  of  new  and  more  accurate  density  fimctionals;  (2)  the  increasing  versatility,  efficiency,  and 
availability  of  DFT  codes;  and  (3)  the  superior  ratio  of  accuracy  to  effort  exhibited  by  DFT  computations 
relative  to  other  ab  initio  methodologies.  It  has  been  clear  for  some  time  that  harmonic  force  fields 
calculated  via  DFT  are  substantially  more  accurate  than  self-consistent  field  (SCF)  force  fields,  and  are 
comparable  in  accuracy  to  force  fields  calculated  using  MoUer-Plesset  second  order  perturbation  theory 
(MP2)  (Ziegler  1991;  Labanowski  and  Andzelm  1991;  Andzelm  and  Wimmer  1992;  Johnson,  GUI,  and 
Pople  1993;  Handy,  Murray,  and  Amos  1993).  Very  recently  there  have  been  two  very  important 
developments  that  greatly  enhance  the  utility  of  DFT  in  predictions  of  vibrational  spectra.  First,  a  new 
class  of  density  functionals  has  been  introduced — the  so-caUed  hybrid  functionals — which  are  significantly 
more  accurate  than  prior  functionals  (Becke  1993).  Second,  efficient  analytical  derivative  techniques  have 
been  implemented  for  the  calculation  of  second  derivatives  of  the  energy  with  respect  to  nuclear 
displacement  (Johnson  and  Frisch  1993,  1994).  As  a  results  of  these  developments,  it  is  now  possible  to 
calculate  harmonic  force  fields  of  an  accuracy  comparable  to  that  of  MP2  force  fields  with  substantially 
less  computational  effort.  In  this  report,  we  document  the  practicabUity  and  illustrate  the  accuracy  of 
calculations  of  VCD  spectra  based  on  DFT  force  fields  calculated  using  a  hybrid  density  functional.  The 
molecules,  4-methyl-2-oxetanone,  i,  6,8-dioxabicyclo[3.2.1]  octane,  2,  and  l,7,7-trimethylbicyclo[2.2.1] 
heptan-2-one  (camphor)k,  3,  are  shown  in  Figure  1. 


Figure  1.  Molecules  studied  in  this  woric;  1,  4-methvl-2-oxetanone;  2,  6.8-dioxabicvclor3.2.11  octane; 
3.  and  1.7.7-trimethvlbicvclor2.2.11  heptan-2-one. 
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Why  is  the  advance  we  report  so  important  for  VCD  spectroscopy?  Since  the  discoveiy  of  VCD  in 
the  1970s  (Holzwarth  et  al.  1974;  Nafie,  Cheng,  and  Stephens  1975;  Nafie,  Keiderling,  and  Stephens  1976; 
Stephens  and  Qaric  1979),  the  primary  justification  for  continued  study  of  VCD  has  been  its  potential 
utility  in  the  elucidation  of  the  stereochemistry  of  chiral  molecules,  both  natural  and  synthetic.  However, 
to  date,  VCD  spectroscopy  has  made  very  little  impact.  While  the  measurement  of  VCD  spectra  has 
become  almost  routine  (Nafie  1988;  Keiderling  1990),  the  extraction  of  structural  information  from  spectra 
has  not.  The  latter  requires  a  methodology  for  accurately  predicting  VCD  spectra.  This,  in  turn,  requires 
ab  initio  computational  methods.  While  much  progress  has  been  made  since  the  earliest  ab  initio 
calculations  of  VCD  spectra  (Lowe,  Stephens,  and  Segal  1986;  Jalkanen  et  al.  1987, 1988;  Kawiecki  et  al. 

1988) ,  the  calculation  of  accurate  harmonic  force  fields  has  remained  problematical.  It  was  clear  a  long 
time  ago  that  SCF  force  fields  are  insufficiently  accurate  to  provide  useful  predictions  of  VCD  spectra 
(except  in  the  case  of  a  few,  very  small  symmetrical  molecules).  MP2  force  fields  provide  much  more 
accurate  spectra  (Amos,  Handy,  and  Palmieri  1990;  Stephens  et  al.  1993;  Stephens  et  al.  1994).  However, 
the  computational  demands  of  MP2  force  fields  increase  very  rapidly  with  increasing  molecular  size,  and 
useful  calculations  on  such  molecules  as  camphor  (Devlin  and  Stephens  1994)  have  only  been  practicable 
on  very  large  supercomputers.  As  we  demonstrate  here,  DFT  provides  a  way  out  of  this  dilemma.  It 
includes  correlation  but  without  the  extra  cost  of  MP2  calculations.  As  a  result,  the  way  is  now  clear  to 
the  routine  calculation  of  VCD  spectra  with  useful  accuracy  on  a  wide  range  of  molecules  up  to  a  very 
considerable  size.  Widespread  application  of  VCD  spectroscopy  is  now,  finally,  close  at  hand. 

2.  METHODS 

Dipole  and  rotational  strengths,  Dj  and  Rj,  of  a  fundamental  transition  in  the  ith  normal  mode  of 
energy  h©j  are  calculated  within  the  harmonic  approximation  using  the  equations  (Stephens  1985,  1987, 

1989) : 
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Normal  coordinates  Qj  are  given  by 
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where  are  Cartesian  displacement  coordinates.  Atomic  polar  tensors  (APTs)  are  given  by 
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where  'Fq  is  the  ground  electronic  state  wavefunction.  Atomic  axial  tensors  (AATs)  are  given  in  the 
distributed  origin  (DO)  gauge  (Stephens  1987;  Jalkanen  et  al.  1988;  Stephens  et  al.  1989, 1990)  and  with 
respect  to  molecular  origin  O  by 
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and  is  calculated  with  the  origin  at  the  equilibrium  position  of  nucleus  k,  r“.  The  "local"  AATs, 

are,  in  turn,  obtained  via  the  identity 
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where  Ej^p(7c)  is  the  electronic  component  of  the  APT  in  the  momentum  representation  (Stephens  et  al. 
1990;  Amos,  Jalkanen,  and  Stephens  1988).  The  DO  gauge  provides  origin-independent  rotational 
strengths. 

For  the  present  work,  DFT  Cartesian  harmonic  force  fields  have  been  calculated  for  1-3  using 
analytical  derivative  techniques  via  the  program  GAUSSIAN  92/DFT  (Frisch  et  al.  1992,  1993).  APTs 
are  calculated  simultaneously.  Three  density  functionals  have  been  employed;  LSDA,  BLYP,  and  the 
hybrid  functional  Becke3LYP  (henceforth  B3LYP).  The  LSDA  functional  (also  referred  to  elsewhere  as 
SVWN)  uses  the  standard  local  exchange  functional  (Becke  1989)  and  the  local  correlation  functional  of 
Vosko,  Wilk,  and  Nusair  (1980).  The  BLYP  functional  combines  the  standard  local  exchange  functional 
with  the  gradient  correction  of  Becke  (1989)  and  uses  the  Lee- Yang-Parr  (1988)  correlation  functional 
(which  also  includes  density  gradient  terms).  The  B3LYP  functional  is  a  hybrid  of  exact  (Hartree-Fock) 
exchange  with  local  and  gradient-corrected  exchange  and  correlation  terms,  as  first  suggested  by  Becke 
(1993a).  The  exchange-correlation  functional  proposed  and  tested  by  Becke  (1993b)  was 
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Here  AE^®®^  is  Becke’s  gradient  correction  to  the  exchange  functional  and  AE/^^^  is  the  Perdew-Wang 
gradient  correction  to  the  correlation  functional  (Perdew  1991).  Becke  (1993b)  suggested  coefficients 
a^  =  0.2,  a^  =  0.72  and  a^  =  0.81  based  on  fitting  to  heats  of  formation  of  small  molecules.  Only  single¬ 
point  energies  were  involved  in  the  fit;  no  molecular  geometries  or  frequencies  were  used.  The  B3LYP 
functional  in  GAUSSIAN  92/DFT  uses  the  values  of  a^,  a^,  and  a^.  suggested  by  Becke,  but  uses  LYP  for 
the  correlation  functional.  Since  LYP  does  not  have  an  easily  separable  local  component,  the  VWN  local 
correlation  expression  has  been  used  to  provide  the  different  coefficients  of  local  and  gradient  corrected 
correlation  functionals: 
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The  standard  fine  grid  in  GAUSSIAN  92/DFT  (Trucks  and  Frisch,  to  be  published)  was  used  in  all 
DFT  calculations.  This  grid  was  produced  from  a  basic  grid  having  75  radial  shells  and  302  angular 
points  per  radial  shell  for  each  atom,  and  reducing  the  number  of  angular  points  for  different  ranges  of 
radial  shells,  leaving  about  7,000  points  per  atom  while  retaining  similar  accuracy  to  the  original 
(75,302)  grid.  Becke’s  (1988)  numerical  integration  techniques  were  employed.  In  addition,  in  the  case 
of  i,  Cartesian  harmonic  force  fields  and  APTs  were  also  calculated  at  the  SCF  and  MP2  levels  of 
approximation  using  analytical  derivative  techniques  via  GAUSSIAN  92  (Frisch  et  al.  1992, 1993).  Irocal 
AATs  were  calculated  for  1-3  at  the  SCF  level  of  approximation  using  analytical  derivative  techniques 
(Amos,  Jalkanen,  and  Stephens  1988;  Amos  et  al.  1987)  via  CADPAC  (Amos  1993).  DO  gauge  AATs 
using  SCF,  MP2,  and  DFT  APTs  are  referred  to  as  SCF,  semi-MP2,  and  semi-DFT  AATS,  respectively. 
Note  that  the  use  of  the  DO  gauge  permits  the  partial  inclusion  of  correlation  in  AATs  when  correlated 
APTs  are  available. 

Frequencies,  dipole  strengths  and  rotational  strengths  derived  from  harmonic  force  fields,  APTs  and 
AATs  were  used  to  synthesize  absorption  and  VCD  spectra  using  Lorentzian  band  shapes  (Kawiecki 
1988). 

Calculations  used  3-21G  (Hehre  et  al.  1986),  6-31G*  (Hehre  et  al.  1986),  and  [5s4p2d/3s2p]  TZ/2P 
(Stephens  et  al.  1990)  basis  sets. 

3.  RESULTS 

We  begin  by  comparing  xmpolarized  absorption  and  VCD  spectra  calculated  using  6-31G*  B3LYP 
DFT  force  fields  to  experimental  spectra.  The  6-31G*  basis  set  (Hehre  et  al.  1986)  is  widely  used  in  ab 
initio  computational  studies  of  organic  molecules,  especially  in  calculations  including  correlation  (e.g., 
MP2).  We  examine  only  the  mid-IR  spectral  region  and  exclude  the  C-H  stretching  region  from  the 
discussion.  Without  explicit  consideration  of  an  harmonicity,  calculations  of  C-H  stretching  spectra  are 

futile.  Experimental  mid-IR  VCD  spectra  are  limited  to  frequencies  >  650  cm“^  (the  current  lower  limit 

of  existing  instrumentation).  We  limit  our  discussion  of  unpolarized  absorption  spectra  to  the  region  over 
which  VCD  spectra  have  been  measured. 
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Figure  2a-d  shows  calculated  and  experimental  mid-IR  spectra  for  1.  The  lower  frequency  limit  of 
the  experimental  spectra  is  690  cm~^  Comparison  of  calculated  and  experimental  absorption  spectra 
yields  the  unambiguous  assignment  of  fundamentals  7-23  shown  in  Figures  2a  and  2b.  The  assignment 
is  further  detailed  in  Table  1.  Calculated  frequencies  are  in  all  cases  greater  than  experimental 
frequencies;  the  average  deviation  of  calculated  and  experimental  frequencies  is  39  cm“^  (3.3%)  and  the 
maximum  deviation  is  74  cm~^  (5.1%). 

Calculated  and  experimental  VCD  spectra  are  compared  in  Figures  2c  and  2d.  VCD  is  observed  for 
all  fundamentals  excepting  22  and  23.  The  agreement  of  calculated  and  observed  VCD  intensities  is  quite 
good,  being  worst  for  fundamentals  12  and  18.  The  sign  is  correctly  predicted  for  aU  bands  except  18. 

Figure  3  shows  calculated  and  experimental  mid-IR  spectra  for  2.  The  lower  frequency  limit  of  the 
experimental  spectra  is  800  cm~^.  Comparison  of  calculated  and  experimental  absorption  spectra  yields 
the  assignment  of  the  fundamentals  10-38  shown  in  Figures  3a  and  3b.  The  assignment  is  unambiguous 
for  the  fundamentals  in  the  range  800-1,300  cm~^;  above  1,300  cm"\  the  lower  intensity  and  greater 
congestion  of  the  spectrum,  together  with  the  greater  probability  of  contributions  from  nonfundamentals, 
lessens  the  confidence  level  somewhat.  The  assigiunent  is  further  detailed  in  Table  2.  Calculated 
frequencies  are  in  aU  cases  greater  than  experimental  frequencies;  the  average  deviation  of  calculated  and 
experimental  frequencies  is  36  cm"^  (2.9%)  and  the  maximum  deviation  is  75  cm"^  (5.2%). 

Comparison  of  the  calculated  and  experimental  VCD  spectra  is  complicated  by  the  much  lower  signal- 
to-noise  ratio  of  the  experimental  VCD  spectrum,  in  comparison  to  the  absorption  spectrum. 
Fimdamentals  12-15,  17-20,  23,  24,  26,  28-33,  and  38  exhibit  clearly  defined  VCD.  For  these  bands, 
the  calculated  VCD  agrees  in  sign  with  the  exception  of  28  and  29.  Below  1,300  cm“^  the  agreement 
of  calculated  and  observed  VCD  intensities  is  quite  good;  the  agreement  is  worst  for  fundamental  26.  The 
calculation  significantly  underestimates  the  intensities  of  the  bands  in  the  region  1,300-1,400  cm“^ 

Figure  4  shows  calculated  and  experimental  mid-IR  spectra  for  3.  The  lower  fiiequency  limit  of  the 
experimental  spectra  is  900  cm“^  Comparison  of  calculated  and  experimental  absorption  spectra  yields 
the  assignment  of  the  fundamentals  24-58  shown  in  Figures  4a  and  4b.  Fundamentals  29-40,  42-44, 
46-50  are  clearly  resolved  and  unambiguously  assignable.  Fundamentals  41  and  45  are  weak  and  not 
clearly  resolved.  The  regions  900-950  cm“^  and  1,440-1,500  cm"^  are  more  congested  and  the 
assignments  of  fundamentals  24-28  and  51-58  are  less  certain.  The  assignment  is  further  detailed  in 
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Shanes  with  v  =  4.0  cm~^  for  all  bands.  Note  that  some  of  the  differences  between  calculated 


Table  1.  4-Methyl-2-Oxetanone 
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and  d)  use  6-31G*  B3LYP/DFT  force  field  and  APTs;  AATs  are  6-31G*  semi-B3LYP/DFT. 


’=  4.0  cm  ^  for  all  bands.  The  assignment  of  fundamentals  is  indicated;  see  also  Table  2. 


Table  2.  6,8-Dioxabicyclo[3.2.1]Octane® 


Calculation‘s 

Experiment‘d 

Fundamental^ 

V 

D 

R 

V 

D 

R 

38 

1,555 

4.1 

7.2 

1,484 

13.1 

19.4 

37 

1,532 

12.0 

1,458 

20.7 

36 

1,513 

5.4 

1.0 

1,438 

11.9 

35 

1,506 

4.8 

-0.4 

1,435 

5.8 

1,430 

6.7 

34 

1,415 

17.2 

-0.3 

1,365 

11.1 

33 

1,412 

25.0 

-6.5 

1,362 

14.9 

-18.8 

32 

1,401 

14.6 

2.4 

1,343 

21.8 

24.9 

31 

1,386 

11.8 

-6.3 

1,338 

11.4 

-29.8 

30 

1,376 

39.4 

7.0 

1,331 

37.6 

52.3 

29 

1,359 

0.8 

0.7 

1,315 

5.8 

-15.2 

28 

1,349 

30.9 

-0.6 

1,310 

28.1 

11.6 

27 

1,314 

1.7 

0.2 

1,276 

0.8 

26 

1,278 

12.8 

-2.2 

1,240 

11.5 

-16.4 

25 

1,243 

4.2 

0.2 

1,211 

24 

1,214 

25.9 

5.7 

1,182 

34.4 

23 

1,187 

112.6 

8.6 

1,157 

126.3 

13.6 

22 

1,149 

328.7 

-21.7 

1,124 

171.5 

1,119 

205.9 

21 

1,120 

20.4 

2.5 

1,089 

33.4 

-9.6 

20 

1,095 

18.6 

18.0 

1,076 

17.0 

19.3 

19 

1,053 

78.0 

30.0 

1,033 

77.4 

20.3 

18 

1,047 

149.9 

-36.0 

1,022 

190.8 

-62.3 

17 

1,023 

180.7 

-29.0 

993 

248.5 

-80.0 

16 

979 

25.0 

-7.3 

962 

32.6 

6.4 

15 

953 

117.3 

71.5 

939 

120.2 

88.0 

14 

911 

64.5 

20.5 

893 

120.1 

63.3 

13 

903 

103.2 

-27.9 

882 

216.0 

-74.0 

12 

878 

145.3 

29.4 

857 

174.6 

32.5 

11 

842 

8.9 

-2.3 

832 

26.1 

10 

815 

16.7 

-4.5 

809 

21.2 

®  V  in  cm  D  in  10“^®  esu^  cm^,  R  in  10~^  esu^  cm^.  Rotational  strengths  are  for  the  (IR,  5S)  -  (+)  enantiomer. 
Fundamentals  are  numbered  in  order  starting  from  the  lowest  frequency.  Nine  fundamentals  of  2  lie  below  800  cm”*. 
6-31G*  B3LYP/DFT  harmonic  force  field  and  APTs;  6-31G*  semi-B3LYP/DFT  AATs. 

^  Eggimann,  Shaw,  and  Wieser  (1991)  and  Eggimarm  et  al.  (1993).  Note  that  our  assignment  of  fundamentals  differs  from  that 
of  Wieser  and  coworkers.  Where  oiu  assignments  are  not  unambiguous,  all  alternative  possibilities  are  indicated. 
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fundamentals  is  indicated:  see  also  Table  3. 
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Table  3.  Calculated  frequencies  are  in  all  cases  greater  than  experimental  frequencies.  The  average 
deviation  of  calculated  and  experimental  frequencies  is  41  cm“^  (3.2%);  the  maximum  deviation  is 
79  cm"^  (5.3%). 

Comparison  of  the  calculated  and  experimental  VCD  spectra  shown  in  Figures  4c  and  4d  is  again 
complicated  by  the  lower  signal-to-noise  ratio  of  the  experimental  VCD  spectrum.  Fundamentals  24,  25, 
27,  30-32,  36,  37,  39, 42, 44,  46-50,  and  57  exhibit  clearly  defined  VCD.  For  these  bands,  the  calculated 
VCD  agrees  in  sign  in  every  case  with  the  exception  of  fundamental  30.  The  agreement  of  calculated  and 
observed  VCD  intensities  is  quite  good,  being  worst  for  fundamentals  32, 44,  and  50.  Despite  the  greater 
complexity  of  the  regions  900-950  cm“^  and  1,440-1,500  cm"^  calculated  and  experimental  VCD  spectra 
are  in  quite  good  agreement. 

The  level  of  agreement  of  calculated  and  experimental  mid-IR  absorption  and  VCD  spectra  exhibited 
in  Figures  2-4  is  impressive.  However,  the  agreement  is  not  perfect.  The  possible  sources  of  error  are: 
(1)  errors  in  the  force  field,  APTs,  and  AATs  due  to  the  incompleteness  of  the  6-31G*  basis  set;  (2)  errors 
in  the  local  AATs  due  to  the  absence  of  correlation;  (3)  errors  in  the  force  field  and  APTs  due  to  the 
inexactness  of  the  B3LYP  density  functional;  and  (4)  the  absence  of  an  harmonicity  and  condensed 
phase/solvent  perturbations  in  the  calculational  formalism. 

Basis  set  error  can  be  determined  from  calculations  using  large  basis  sets.  In  the  case  of  j[,  we  have 
carried  out  calculations  using  a  [5s4p2d/3s2p],  TZ/2P,  basis  set,  with  the  results  given  in  Figure  5  and 
Table  1.  The  agreement  of  calculated  and  experimental  frequencies  is  improved:  the  average  and 
maximum  deviations  are  now  25  cm“^  (2.0%)  and  48  cm~^  (3.3%),  respectively.  Both  calculated 
absorption  and  VCD  intensities  are  in  better  agreement  with  experiment.  The  improvement  is  more 
marked  in  the  case  of  the  VCD  spectrum:  fundamental  18  is  now  correct  in  sign  and  the  intensities  of 
fundamentals  1 1  and  12  are  in  significantly  improved  agreement  with  experiment.  The  greater  basis  set 
dependence  of  the  calculated  VCD  spectra  than  of  the  calculated  absorption  spectra  is  attributable  to  the 
greater  basis  set  error  in  the  local  AATs  than  in  the  force  field  and  APTs.  This  is  verified  by  calculations 
using  the  6-3 IG*  force  field  and  APTs  together  with  TZ/2P  local  AATs,  the  results  of  which  are  also 
given  in  Figure  5. 

We  are  not  yet  able  to  examine  the  magnitudes  of  errors  originating  in  the  absence  of  correlation  in 
local  AATs  by  canying  out  calculations  in  which  correlation  is  included.  However,  we  can  examine  the 
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Table  3.  Camphoi^ 


^  V  in  cm"’,  D  in  lO"^  esu^  cm^,  R  in  10"^  esu^  cm’’;  anisotropy  ratio  g  in  10“^.  Rotational  strengths  are  for  the 
(IR,  4R)  -  (+)  enantiomer. 

^  Fundamentals  are  numbered  in  order  starting  from  the  lowest  frequency.  Twenty-three  fundamentals  of  3  lie  below 
900  cm"’. 

^  6-31G*  B3LYP/DFT  harmonic  force  field  and  APTs;  6-31G*  semi-B3LYP/DFT  AATs. 

^  Devlin  and  Stephens  (1994)  and  Nafie  (1984). 
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Lorentzian  fitting  of  the  experimental  spectra  (Figures  2a  and  2c)  followed  bv  elimination  of 


all  nonfundamental  bands.  Calculated  spectra  (b,  c.  e.  f,  and  g)  use  B3LYP/DFT  force  field 


and  APTs;  AATs  are  semi-B3LYP/DFT.  In  b  and  e,  the  basis  set  is  TZ/2P:  in  c  and 


6-31G*.  In  f.  the  force  field  and  APTs  are  6-31 G*  while  the  local  AATs  are  T 


4.0  cm  ^  for  aU  bands  in  aU  calculated  soectra. 
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results  of  calculations  in  which  correlation  is  not  included  in  the  APTs.  It  is  not  unreasonable  to  expect 
that  the  inclusion  of  correlation  in  APTs  and  in  AATs  wiU  give  rise  to  contributions  to  rotational  strengths 
of  comparable  order  of  magnitude.  The  results  for  1  obtained  using  6-31G*  SCF  APTs  in  place  of  631G* 
B3LYP/DFT  APTs  are  given  in  Figure  6  and  Table  1.  Both  absorption  and  VCD  spectra  are  qualitatively 
unaffected.  Quantitatively,  the  changes  are  larger  in  the  absorption  intensities,  which  are  generally 
increased  by  the  substitution  of  SCF  for  B3LYP/DFT  APTs.  It  is  reasonable  to  conclude  that  inclusion 
of  correlation  in  the  calculation  of  local  AATs  will  not  cause  major  qualitative  changes  in  the  calculated 
VCD  spectrum. 

At  this  time,  there  does  not  appear  to  be  any  density  functional  available  whose  accuracy  is 
significanfly  greater  than  that  of  the  B3LYP  functional.  We  cannot  therefore  directly  examine  the  errors 
in  B3LYP/DFT  calculations  by  comparison  to  more  accurate  calculations.  However,  we  can  examine  the 
sensitivity  of  the  calculated  spectra  to  alternative,  less  accurate  choices  of  density  functional.  The  results 
for  1  obtained  using  two  widely  employed  functionals — LSDA  and  BLYP — are  compared  to  those 
obtained  using  B3LYP  in  Figure  7.  Both  functionals  give  absorption  and  VCD  spectra  substantially 
different  from  and  in  much  worse  agreement  with  experiment  than  given  by  the  B3LYP  functional.  It  is 
clear  that  the  LSDA  and  BLYP  harmonic  force  fields  are  significantly  less  reliable  than  the  B3LYP  force 
field.  In  turn,  it  is  certainly  possible  that  some  of  the  errors  in  predicted  B3LYP  absorption  and  VCD 
spectra  are  attributable  to  residual  inexactitude  in  the  B3LYP  functional. 

The  results  shown  in  Figure  5  and  Table  1  demonstrate  that  the  6-3 IG*  B3LYP  force  field  is  not 
substantially  improved  if  a  much  larger  basis  set  is  substituted  for  6-31G*.  It  is  therefore  of  interest  to 
examine  the  consequences  of  using  a  smaller  basis  set.  Spectra  calculated  for  !_  using  the  3-21G  basis  set 
are  shown  in  Figure  8.  Both  absorption  and  VCD  spectra  are  substantially  different  from  those  obtained 
at  the  6-3 IG*  basis  set  level  and  in  much  worse  agreement  with  experiment.  It  is  clear  that  a  basis  set 
at  least  as  large  as  6-31G*  is  mandatoiy  if  useful  agreement  with  experiment  is  to  be  obtained. 

Finally,  in  Figure  9,  we  compare  the  spectra  for X  calculated  using  the  6-31G*  B3LYP/DFT  force  field 
to  the  spectra  obtained  using  6-31G*  SCF  and  6-31G*  MP2  force  fields.  APTs  are  calculated  at  the  same 
level  as  the  force  field;  local  AATs  are  in  all  cases  calculated  at  the  6-3 IG*  SCF  level.  The  MP2 
methodology  yields  spectra  quite  similar  to  those  obtained  using  DFT  with  the  B3LYP  functional  and  in 
slightly  worse  agreement  with  experiment.  The  most  noticeable  difference  is  in  the  relative  absorption 
intensities  of  fundamentals  14  and  15.  In  contrast,  the  SCF  methodology  yields  quite  different  spectra  in 
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Figure  6.  Mid-IR  absorption  and  VCD  spectra  of  R  -  (+)  -  1.  Spectra  a  and  d  are  the  fundamental 
experimental  spectra  from  Figures  5a  and  5d.  Calculated  spectra  fb.  c,  e.  and  f)  use  the  6-31G* 
B3LYP/DFT  force  field.  In  b  and  e.  the  APTs  are  6-3 IG*  B3LYP/DFT  and  the  AATs  ate 
6-31G*  semi-B3LYP/DFT.  In  c  and  f,  the  APTs  and  AATs  are  6-31G*  SCF.  y=  4.0  cm~^ 
for  all  bands  in  all  calculated  spectra. 
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Figure  8.  Mid-IR  absorption  and  VCD  spectra  of  R  -  (+)  -  1.  Spectra  a  and  d  are  the  fundamental 
experimental  spectra  from  Figures  5a  and  5d.  Calculated  spectra  (b.  c,  e,  and  f)  use 
B3LYP/DFT  force  field  and  APTs;  AATs  are  semi-B3LYP/DFT.  In  b  and  e.  the  basis  set  is 
6-31G*:  in  c  and  f,  it  is  3-21G.  Y  =  4.0  cm~^  for  all  bands  in  all  calculated  spectra. 


18 


wavenumbers 


very  much  worse  agreement  with  experiment  The  importance  of  including  correlation  in  the  calculation 
of  harmonic  force  fields  is  forcefully  demonstrated  by  these  results.  At  the  same  time,  the  B3LYP/DFT 
method  is  clearly  shown  to  be  highly  competitive  with  the  MP2  methodology  in  including  correlation. 

Analogous  calculations  for  2  and  3  to  those  for  X  presented  in  Figures  6-9  confirm  the  generality  of 
the  conclusions  arrived  at  above  in  the  case  of  1.  These  calculations  will  be  reported  elsewhere.  (Note 
that  6-31G*  MP2  calculations  have  already  been  published  for  3  [Devlin  and  Stephens  1994]). 

In  the  case  of  2,  a  second  conformation  is  possible  in  which  the  six-membered  ring  is  in  a  boat,  rather 
than  a  chair,  conformation.  Ab  initio  calculations  (6-31G*  SCF,  MP2,  and  B3LYP/DFT)  predict  a 
4-6  kcalAnole  energy  difference.  Mid-IR  absorption  and  VCD  spectra  predicted  for  the  "boat" 
configuration  of  2  are  compared  to  the  experimental  spectra  in  Figure  10.  In  contrast  to  the  spectra 
predicted  for  the  "chair"  conformation,  the  agreement  with  the  experimental  spectra  is  very  poor.  The 
difference  is  especially  dramatic  in  the  case  of  the  VCD  spectra.  As  expected,  the  "chair"  conformation 
of  2  is  unquestionably  that  existing  experimentally. 

4.  DISCUSSION 

We  have  demonstrated  that  calculations  of  the  mid-IR  absorption  and  VCD  spectra  of  X-3  based  on 
6-31G*  B3LYP/DFT  harmonic  force  fields  are  (1)  in  impressive  agreement  with  experiment; 

(2)  substantially  more  accurate  than  calculations  using  the  LSDA  and  BLYP  density  functionals; 

(3)  comparable  in  accuracy  to  calculations  based  on  6-3 IG*  MP2  force  fields,  while  much  less 
computationally  demanding.  Substitution  of  TZ/2P  and  3-2 IG  basis  sets  for  6-3 IG*  leads  to  spectra  of 
very  similar  and  much  worse  accuracy,  respectively. 

While  the  considerable  accuracy  of  DFT  calculations  using  "modem”  density  functionals  has  been 
widely  reported  (Ziegler  1991;  Labanowski  and  Andzelm  1991;  Andzelm  and  Wimmer  1992;  Johnson, 
GiU,  and  Pople  1993;  Handy,  Murray,  and  Amos  1993),  the  sensitivity  of  predicted  vibrational  spectra  to 
the  choice  of  density  functional  was  not  clear  prior  to  this  work.  Comparisons  of  harmonic  frequencies 
calculated  using  various  functionals,  including  LSDA  (=SVWN)  and  BLYP,  to  experimental  harmonic 
frequencies  did  not  suggest  large  differences  in  accuracy  (Johnson,  Gill,  and  Pople  1993).  However, 
dipole  and  rotational  strengths  depend  on  vibrational  coordinates,  which  are  much  more  sensitive  to  the 
accuracy  of  the  force  field  than  are  vibrational  frequencies.  It  is  therefore  not  surprising  to  find  that 
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calculated  spectra  are,  in  fact,  quite  sensitive  to  the  choice  of  density  functional.  The  really  new  result 
is  that  the  hybrid,  B3LYP,  functional  gives  spectra  of  impressive  accuracy  while  the  LSDA  and  BLYP 
functionals  do  not.  Our  work  thus  defines  a  threshold  level  of  accuracy  for  density  functionals  in  future 
calculations  of  mid-IR  vibrational  spectra.  At  the  same  time,  it  is  clear  that  further  improvements  in 
density  functionals  may  well  lead  to  significantly  improved  vibrational  spectra.  B3LYP  could  be  replaced 
as  the  functional  of  choice  in  the  near  future! 

Our  woik  also  further  illuminates  the  relative  accuracies  of  DFT  and  MP2  harmonic  force  fields.  At 
the  B3LYP  level,  DFT  and  MP2  force  fields  are  very  similar  in  accuracy  when  this  is  gauged  by  the 
comparison  of  calculated  and  experimental  mid-IR  spectra  for  1-3.  The  LSDA  and  BLYP  DFT  force 
fields  are  clearly  inferior  to  MP2  force  fields.  This  latter  finding  is  surprising:  it  has  been  reported  that 
such  functionals  yield  harmonic  frequencies  comparable  in  accuracy  with  MP2  calculations  (Johnson,  Gill, 
and  Pople  1993).  It  is  possible  that  this  difference  reflects  the  focus  of  our  work  on  mid-IR  spectra,  to 
the  exclusion  of  hydrogenic  stretching  and  low-frequency  spectral  regions,  while  comparisons  of  calculated 
and  experimental  frequencies  for  small  molecules  (the  only  ones  for  which  harmonic  frequencies  are 
known)  include  all  frequencies  and  are  heavily  weighted  by  hydrogenic  stretching  frequencies. 

The  conclusions  we  have  arrived  at  have  made  use  of  both  absorption  and  VCD  spectra  in  comparing 
theory  and  experiment  While  the  same  conclusions  would  have  been  reached  if  absorption  or  VCD 
spectra  had  been  used  alone,  the  combined  use  of  both  spectra  significantly  enhances  the  reliability  of  the 
assignments  of  fundamentals  and  of  the  evaluation  of  the  relative  accuracies  of  different  calculations.  The 
usefulness  of  VCD  spectra  for  this  purpose  is  of  course  dependent  on  the  reliability  with  which  VCD 
intensities  can  be  predicted.  This  is  lower  than  the  reliability  of  predicted  absorption  intensities  since 

(1)  absorption  intensities  require  only  APTs,  while  VCD  intensities  require  both  APTs  and  AATs,  and 

(2)  AATs  are  less  accurately  calculated  due  to  (a)  the  absence  of  correlation  in  the  calculation  of  local 
AATs  and  (b)  the  greater  basis  set  error  in  local  AATs  than  in  APTs.  As  expected,  the  agreement  of 
calculated  and  experimental  VCD  spectra  is  generally  less  perfect  than  in  the  case  of  absorption  spectra. 
Nevertheless,  it  is  clear  that  the  6-31G*  semi-B3LYP/DFT  AATs  used  in  this  work  are  of  sufficient 
reliability  to  successfully  predict  the  qualitative  pattern  of  VCD  intensities  with  a  very  small  number  of 
exceptions,  and  to  support  the  interrogation  of  alternative  force  fields  in  a  meaningful  way. 

Currently,  the  practical  limits  in  calculations  of  harmonic  force  fields  are  imposed  by  the  disk  storage 
requirements.  Analytical  derivative  calculations  at  the  MP2  level  require  much  more  disk  space  than  do 
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SCF  calculations.  DFT  calculations  have  identical  requirements  to  SCF  calculations.  SCF  and  DFT 
calculations  can  therefore  be  carried  out  for  much  larger  molecules  than  are  accessible  to  MP2 
calculations.  Currently,  with  the  machines  available  to  us  (which  include  a  CRAY  C90)  MP2  6-31G* 
calculations  on  3  are  close  to  the  limit  of  practicability.  In  contrast,  DFT  6-3 IG*  calculations  wiU  be 
straightforward  for  molecules  much  larger  than  3.  Our  results  here  lead  to  the  conclusion  that  these 
calculations  wiU  be  of  useful  accuracy  as  long  as  the  B3LYP  (or  better)  density  functional  is  used. 

5.  CONCLUSION 

The  implications  for  the  utilization  of  VCD  spectroscopy  are  obvious  and  dramatic.  The  reliable 
prediction  of  VCD  spectra  wiU  be  routine  for  molecules  much  larger  than  3.  In  turn,  the  elucidation  of 
molecular  stereochemistry  from  VCD  spectra  wiU  be  straightforward.  AU  that  remains  is  to  choose 
interesting  and  important  problems  to  study. 
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DIVISION  8354 
ATTN  S  JOHNSTON 
PMATTERN 

D  STEPHENSON 

SANDIA  NATIONAL  LABORATORIES 
LIVERMORE  CA  94550 

1  BRIGHAM  YOUNG  UNIVERSITY 

ATTN  M  W  BECKSTEAD 
DEPT  OF  CHEMICAL  ENGINEERING 
PROVO  UT  84058 

1  CALIFORNIA  INSTITUTE  OF  TECH 

JET  PROPULSION  LABORATORY 
ATTN  L  STRAND  MS  125  224 
4800  OAK  GROVE  DRIVE 
PASADENA  CA  91109 

1  CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 

ATTN  F  E  C  CULICK  MC  301  46 
204KARMANLAB 
PASADENA  CA  91125 

1  UNIVERSITY  OF  CALIFORNIA 

LOS  ALAMOS  SCIENTIFIC  LAB 
P  O  BOX  1663  MAIL  STOP  B216 
LOS  ALAMOS  NM  87545 

1  UNIVERSITY  OF  CALIFORNIA  BERKELEY 

CHEMISTRY  DEPARMENT 
ATTN  C  BRADLEY  MOORE 
211  LEWIS  HALL 
BERKELEY  CA  94720 

1  UNIVERSITY  OF  CALIFORNIA  SAN  DIEGO 

ATTN  F  A  WILLIAMS 
AMES  BOlO 
LA  JOLLA  CA  92093 
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2  UNIV  OF  CALIFORNIA  SANTA  BARBARA 
QUANTUM  INSTITUTE 

ATTN  K  SCHOFIELD 

M  STEINBERG 

SANTA  BARBARA  CA  93106 

1  UNTV  OF  COLORADO  AT  BOULDER 

ENGINEERING  CENTER 
ATTN  J  DAILY 
CAMPUS  BOX  427 
BOULDER  CO  80309-0427 

3  UNTV  OF  SOUTHERN  CALIFORNIA 
DEPT  OF  CHEMISTRY 

ATTN  R  BEAUDET 
S  BENSON 
c  wnriG 

LOS  ANGELES  CA  90007 

1  CORNELL  UNIVERSITY 

DEPARTMENT  OF  CHEMISTRY 
ATTN  T  A  COOL 
BAKER  LABORATORY 
ITHACA  NY  14853 

1  UNIVERSITY  OF  DELAWARE 

ATTN  T  BRILL 
CHEMISTRY  DEPARTMENT 
NEWARK  DE  19711 

1  UNIVERSITY  OF  FLORIDA 

ATTN  J  WINEFORDNER 
DEPT  OF  CHEMISTRY 
GAINESVILLE  FL  32611 

3  GEORGIA  INSTITUTE  OF  TECHNOLOGY 

SCHOOL  OF  AEROSPACE  ENGINEERING 
ATTN  E  PRICE 
WCSTRAHLE 
BTZINN 

ATLANTA  GA  30332 

1  UNIVERSITY  OF  ILLINOIS 

DEPT  OF  MECH  ENG 
ATTNHKRIER 
144MEB  1206  W  GREEN  ST 
URBANA  IL  61801 

1  THE  JOHNS  HOPKINS  UNIV  CPIA 

ATTN  T  W  CHRISTIAN 
10630  LITTLE  PATUXENT  PKWY  SUITE  202 
COLUMBIA  MD  21044-3200 


1  UNIVERSITY  OF  MICHIGAN 

GAS  DYNAMICS  LAB 
ATTN  G  M  FAETH 
AEROSPACE  ENGINEERING  BLDG 
ANN  ARBOR  MI  48109-2140 

1  UNIVERSITY  OF  MINNESOTA 

DEPT  OF  MECHANICAL  ENGINEERING 
ATTN  E  FLETCHER 
MINNEAPOLIS  MN  55455 

4  PENNSYLVANIA  STATE  UNIVERSITY 
DEPT  OF  MECHANICAL  ENGINEERING 
ATTN K KUO 
MMICCI 
S  THYNELL 
V  YANG 

UNIVERSITY  PARK  PA  16802 

2  PRINCETON  UNIVERSITY 
FORRESTAL  CAMPUS  LIBRARY 
ATTN  K  BREZINSKY 

I  GLASSMAN 
P  O  BOX  710 
PRINCETON  NJ  08540 

1  PURDUE  UNIVERSITY 

SCHOOL  OF  AERONAUTICS  AND  ASTRONAUTICS 

ATTN  J  R  OSBORN 

GRISSOM  HALL 

WEST  LAFAYETTE  IN  47906 

1  PURDUE  UNIVERSITY 
DEPARTMENT  OF  CHEMISTRY 
ATTN  E  GRANT 

WEST  LAFAYETTE  IN  47906 

2  PURDUE  UNIVERSITY 

SCHOOL  OF  MECHANICAL  ENGINEERING 
ATTN  N  M  LAURENDEAU 
S  NB  MURTHY 
TSPC  CHAFFEE  HALL 
WEST  LAFAYETTE  IN  47906 

1  RENSSELAER  POLYTECHNIC  INST 

DEPT  OF  CHEMICAL  ENGINEERING 
ATTN  AFONTUN 
TROY  NY  12181 
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1  STANFORD  UNIVERSITY 

DEPT  OF  MECHANICAL  ENGINEERING 
ATTN  R  HANSON 
STANFORD  CA  94305 

1  UNIVERSITY  OF  TEXAS 

DEPT  OF  CHEMISTRY 
ATTN  W  GARDINER 
AUSTIN  TX  78712 

1  VA  POLYTECH  INST  AND  STATE  UNTV 

ATTN  J  A  SCHETZ 
BLACKSBURG  VA  24061 

1  APPLIED  COMBUSTION  TECHNOLOGY  INC 
ATTN  A  M  VARNEY 

P  O  BOX  607885 
ORLANDO  FL  32860 

2  APPLIED  MECHANICS  REVIEWS 
ASME 

ATTN  R  E  WHITE  &  A  B  WENZEL 
345  E  47TH  STREET 
NEW  YORK  NY  10017 

1  ATLANTIC  RESEARCH  CORP 

ATTN  R  H  W  WAESCHE 
7511  WELLINGTON  ROAD 
GAINESVILLE  VA  22065 

1  TEXTRON  DEFENSE  SYSTEMS 

ATTN  A  PATRICK 
2385  REVERE  BEACH  PARKWAY 
EVERETT  MA  02149-5900 

1  BATTELLE 

TWSTIAC 
505  KING  AVENUE 
COLUMBUS  OH  43201-2693 

1  COHEN  PROFESSIONAL  SERVICES 

ATTN  N  S  COHEN 
141  CHANNING  STREET 
REDLANDS  CA  92373 

1  EXXON  RESEARCH  &  ENG  CO 

ATTN  A  DEAN 
ROUTE  22E 
ANNANDALE  NJ  08801 


1  GENERAL  APPLIED  SCIENCE  LABS  INC 

77  RAYNOR  AVENUE 
RONKONKAMA  NY  11779-6649 

1  GENERAL  ELECTRIC  ORDNANCE  SYSTEMS 

ATTN  JMANDZY 
100  PLASTICS  AVENUE 
PITTSFIELD  MA  01203 

1  GENERAL  MOTORS  RSCH  LABS 
PHYSICAL  CHEMISTRY  DEPARTMENT 
ATTN  T  SLOANE 

WARREN  MI  48090-9055 

2  HERCULES  INC 
ATTN  W  B  WALKUP 
E A  YOUNT 

P  O  BOX  210 

ROCKET  CENTER  WV  26726 

1  HERCULES  INC 

ATTN  R  V  CARTWRIGHT 
100  HOWARD  BLVD 
KENVIL  NJ  07847 

1  ALLIANT  TECHSYSTEMS  INC 

MARINE  SYSTEMS  GROUP 
ATTN  D  E  BRODEN  MS  MN50  2000 
600  2ND  STREET  NE 
HOPKINS  MN  55343 

1  ALLIANT  TECHSYSTEMS  INC 

ATTN  R  E  TOMPKINS 
MN  11  2720 

600  SECOND  ST  NORTH 
HOPKINS  MN  55343 

1  IBM  CORPORATION 

ATTN  A  C  TAM 
RESEARCH  DIVISION 
5600  COTTLE  ROAD 
SAN  JOSE  CA  95193 

1  nr  RESEARCH  INSTITUTE 

ATTN  R  F  REMALY 
10  WEST  35TH  STREET 
CHICAGO  IL  60616 
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1  LOCKHEED  MISSILES  &  SPACE  CO 

ATTN  GEORGE  LO 
3251  HANOVER  STREET 
DEPT  52  35  B204  2 
PALO  ALTO  CA  94304 

1  OUN  ORDNANCE 

ATTN  V  MCDONALD  LIBRARY 
P  O  BOX  222 
ST  MARKS  FL  32355-0222 

1  PAUL  GOUGH  ASSOCIATES  INC 

ATTN  P  S  GOUGH 
1048  SOUTH  STREET 
PORTSMOUTH  NH  03801-5423 

1  HUGHES  AIRCRAFT  COMPANY 

ATTN  T  E  WARD 
8433  FALLBROOK  AVENUE 
CANOGA  PARK  CA  91303 

1  ROCKWELL  INTERNATIONAL  CORP 

ROCKETDYNE  DIVISION 
ATTN  J  E  FLANAGAN  HB02 
6633  CANOGA  AVENUE 
CANOGA  PARK  CA  91304 

1  SCIENCE  APPLICATIONS  INC 

ATTN  R  B  EDELMAN 
23146  CUMORAH  CREST 
WOODLAND  HILLS  CA  91364 

3  SRI  INTERNATIONAL 

ATTN  G  SMITH 
D  CROSLEY 
D  GOLDEN 

333  RAVENSWOOD  AVENUE 
MENLO  PARK  CA  94025 

1  STEVENS  INSTITUTE  OF  TECH 

DAVIDSON  LABORATORY 
ATTN  R  MCALEVY  m 
HOBOKEN  NJ  07030 

1  SVERDRUP  TECHNOLOGY  INC 

LERC  GROUP 

ATTN  R  J  LOCKE  MS  SVR  2 
2001  AEROSPACE  PARKWAY 
BROOK  PARK  OH  44142 


1  SVERDRUP  TECHNOLOGY  INC 

ATTN  J  DEUR 

2001  AEROSPACE  PARKWAY 
BROOK  PARK  OH  44142 

3  THIOKOL  CORPORATION 

ELKTON  DIVISION 
ATTN  R  BIDDLE 
R  WILLER 
TECH  LIB 
P  O  BOX  241 
ELKTON  MD  21921 

3  THIOKOL  CORPORATION 

WASATCH  DIVISION 
ATTN  S  J  BENNETT 
P  O  BOX  524 

BRIGHAM  CITY  UT  84302 

1  UNITED  TECHNOLOGIES  RSCH  CENTER 

ATTN  A  C  ECKBRETH 
EAST  HARTFORD  CT  06108 

1  UNITED  TECHNOLOGIES  CORP 

CHEMICAL  SYSTEMS  DIVISION 
ATTN  R  R  MILLER 
P  O  BOX  49028 
SAN  JOSE  CA  95161-9028 

1  UNIVERSAL  PROPULSION  COMPANY 

ATTN  H  J  MCSPADDEN 
25401  NORTH  CENTRAL  AVENUE 
PHOENIX  AZ  85027-7837 

1  VERTTAY  TECHNOLOGY  INC 

ATTN  E  B  FISHER 
4845  MILLERSPORT  HIGHWAY 
EAST  AMHERST  NY  14051-0305 

1  FREEDMAN  ASSOCIATES 

ATTN  E  FREEDMAN 
2411  DIANA  ROAD 
BALTIMORE  MD  21209-1525 

1  ALLIANT  TECHSYSTEMS 

ATTN  J  BODE 
600  SECOND  ST  NE 
HOPKINS  MN  55343 
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ALLIANT  TECHSYSTEMS 
ATTN  C  CANDLAND 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

ALLIANT  TECHSYSTEMS 
ATTN  L  OSGOOD 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

ALLIANT  TECHSYSTEMS 
ATTN  R  BURETTA 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

ALLIANT  TECHSYSTEMS 
ATTN  R  BECKER 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

ALLIANT  TECHSYSTEMS 
ATTN  M  SWENSON 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

US  ARMY  BENET  LABORATORY 
ATTN  SAM  SOPOK 
SMCAR  CCB  B 
WATERVUET  NY  12189 


ABERDEEN  PROVING  GROUND 


36  DIR  USARL 

ATTN:  AMSRL-WT-P,  A  HORST 
AMSRL-WT-PC, 

R  A  FIFER 
GF  ADAMS 
W  R  ANDERSON 
R  A  BEYER 
S  WBUNTE 
C  F  CHABALOWSKI 
K  P  MC-NEILL  BOONSTOPPEL 
A  COHEN 
RCUMPTON 
R  DANIEL 
DDEVYNCK 
NFFELL 
B  E  FORCH 
JMHEIMERL 
A  JKOTLAR 
MRMANAA 
W  F  MCBRATNEY 
K  L  MCNESBY 
S  VMEDUN 
MS  MILLER 
A  WMIZIOLEK 
S  H  MODIANO 
J  B  MORRIS 
J  E  NEWBERRY 
S  A  NEWTON 
R  A  PESCE-RODRIGUEZ 
BMRICE 
R  C  SAUSA 
M  A  SCHROEDER 
J  A  VANDERHOFF 
M  WENSING 
A  WHREN 
J  M  WIDDER 
C  WILLIAMSON 
AMSRL-CI-CA,  R  PATEL 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-732 _ Date  of  Report  April  192g - 

2.  Date  Report  Received _ _ _ _ _ — 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  i»oject,  or  other  area  of  interest  for  which  die 

report  will  be  used.) _ _ _ _ _ — — 


4.  Specifically,  how  is  flie  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating 
costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. - 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.) _ _ _ _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and 
the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  t^  closed,  and  maiL) 
(DO  NOT  STAPLE) 


DEPARTMENT  OF  THE  ARMY 


OFnCIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  OOOl.APG.MD 


POSTAGE  Wia  BE  PAID  BY  ADDRESSEE 


DIRECTOR 

U.S.  ARMY  RESEARCH  LABORATORY 
ATTN:  AMSRL-WT-PC 

ABERDEEN  PROVING  GROUND,  MD  21005-5066 


